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A control scheme for synchronizing the frequencies of geographically 
separated oscillators connected by communication links consists of averag- 
ing the phases received at each station from remote oscillators, comparing 
the result with the local phase, and applying the filtered error signal as a 
correction to the local oscillator frequency. The system was studied by V. E. 
Benes who found a sufficient condition for the stability of the system using 
advanced mathematical techniques. In this paper, the stability condition 
is derived (for a slightly more general control scheme) using only the transfer 
function concept of linear systems and some properties of determinants. 
A practical difficulty regarding the final frequency of the oscillators is 
discussed and a modification of the control scheme is shown to alleviate the 
difficulty. Also examined are the questions of sensitivity to parameter 
changes, the effect of jitter noise on the performance of the system, and the 
effect of failure of an oscillator or transmission link. 

I. INTRODUCTION 

Consider a network of JV geographically separated stations that are 
connected by directed communication links. A local clock, or oscillator, 
is situated at each station. The problem of synchronizing the frequencies 
of the oscillators is of considerable practical interest for continental 
pulse code modulation (PCM) systems. 

The local oscillators have frequencies which may be altered in pro- 
portion to a control signal. In the absence of external control, each 
oscillator operates at a different frequency. The network is "connected" 
in the sense that from any station to any other station there is either a 
direct transmission link or an indirect path via one or more intermedi- 
ate stations. A fixed time delay is associated with each transmission 
link. 

In an important but unpublished paper, V. E. Benes 1 has examined 
a linear control scheme in which each station receives the phases of 
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neighboring stations, i.e., those stations connected to it by direct trans- 
mission links. The phases are averaged and compared with the local 
phase; the error is filtered and applied as a correction to the frequency 
of the oscillator. Similar schemes were also proposed by Runyon. 2 
Belles' has proved that under suitable conditions the system is stable, 
i.e., the oscillators asymptotically settle to a common frequency and the 
phase differences have finite asymptotic values. He also finds explicit 
formulae for the final frequency and asymptotic phase differences. 
To obtain these results, he resorted to the mathematical techniques of 
renewal theory and Tauberian theory. By assuming the stability of the 
system, as proved by Benes, A. J. Goldstein 3 has rederived the expres- 
sions for final frequency and phase differences in a more direct manner. 
Bonomi, La Marche, and Varaiya 4 improved the treatment of the 
stability problem and suggested some avenues of approach for the 
study of transient response. In each case, the authors relied on the 
mathematical theory of Markoff chains and stochastic matrices. 

M. Karnaugh 5 has formulated a more realistic and more sophisti- 
cated nonlinear control model. Broad stability conditions for this 
model are not yet known; however, certain special cases resemble the 
BeneS model. 

In this paper, the stability conditions and the expression for final 
frequency for a slightly more general version of the Benes model are 
derived in a simple manner using only the transfer function concept of 
linear systems and elementary properties of determinants. This approach 
permits a clearer intuitive understanding and should be readily com- 
prehensible to the non-mathematician. The sensitivity of the system 
to parameter changes is also examined and certain questions regarding 
the final frequency of the oscillators are clarified. 

In Section II we give a formulation of the problem and obtain the 
basic equations describing the system. In Section III certain crucial 
properties of the matrix of averaging coefficients are derived which 
result from the topological constraint that the network is connected. 
Stability is proved in Section IV and an expression for the final fre- 
quency is obtained. Section V considers some practical questions with 
regard to how the final frequency is related to the free-running fre- 
quencies of the oscillators. Section VI examines the questions of sensi- 
tivity to parameter changes, the effect of failure of an oscillator or 
transmission link, and the effect of jitter noise. 

II. FORMULATION 

Let /,- be the frequency of the z'th oscillator in the absence of external 
control, and r,(2) the control signal applied to the ith oscillator at time t . 
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If pi(t) denotes the total cyclical phase of the ith oscillator, then the 
actual frequency at time t is given by 

Pi(t) = ft + n(t) (l) 

where the dot denotes the time derivative. 

The control scheme at the ith oscillator is shown in Fig. 1. The 
phases of all neighboring stations are transmitted to the ith station. 
The transmission delay associated with the path from station j to 
station i is denoted as r# . Each phase received at station i is compared 
with the local phase; the differences are weighted with the nonnegative 
averaging coefficients a,-> and summed. The weighted sum of phase 
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Fig. 1 — Station i of the phase averaging control system. 

differences is applied to the filter with transfer function Hi(s), and the 
filter output is the frequency correction term r,-(<) . Thus, we have 

N 

r t (t) = h>(t) * E a i} {pM - nj) - Pi(t)], (2) 

j=i 

where hi(t) is the impulse response of the filter at the ith station and the 
asterisk denotes convolution. 

We assume the filters have three simple properties: (?) causality, 
i.e., the response at any instant does not depend on the future of the 
input, (it) stability, in the sense that a bounded input always produces 
a bounded output, and (in) positive dc gain, i.e., 

Hi(0) = a, > 0. (3) 

Without loss in generality we may assume that the averaging coeffi- 
cients sum to unity, i.e., 



X an = 1. 



(4) 
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Clearly a scaling of all the coefficients a# for fixed i is equivalent to a 
change in gain factor of the ith filter. 

If there is no direct transmission path from the jth to the ith station, 
then the coefficient a i} - is presumed to be zero. Thus, the N X N ma- 
trix, A , whose ijth element is a,y , contains all the topological informa- 
tion about the network of communication links. In order that mutual 
synchronization be possible, it is certainly necessary that the network be 
connected so that from any station to any other station there is either 
a direct or indirect transmission path. The resulting properties of the 
averaging matrix A imposed by this topological constraint play a 
vital role in the proof of stability for the system. 

In agreement with the Benes model we consider the starting condi- 
tions where the oscillators are assumed to have been free-running 
for an indefinitely long time prior to t = 0, and at t = the control 
signals r, are connected to the oscillators. Thus, we have 

PS) = U + Pi(0), t < (5) 

where Pi(0) is the phase at t = 0, and from (1), (2), and (4), the fre- 
quency of the ith oscillator when the controls are operating is 

Pi(t) = hi(t) * Z atfojit - nj) - pi(t)] + /,-, t ^ 0. (6) 
i 

Equations (5) and (6) for i = 1,2, ■ • • , N completely characterize 
the behavior of the system under the particular starting conditions of 
interest. Taking the ordinary Laplace transform of (6), we obtain 

sPi = Hi E dijPj - HiPi + -fi + p,(0) 4- Qi, (7) 

where 

da = an exp (-stu), 

Qi(s) = HM Z da / Pi (t) exp (-st) dt, 

i J -rii 

and Pi(s) is the Laplace transform of pt(t). The term Q t (s) is the con- 
tribution to the ith oscillator frequency after t = due to the contents 
of the transmission links at t = 0. Using (5), Qi(s) can be evaluated 
explicitly, but for our purposes it is sufficient to note that 

sQi(s) -> as s -> 0. (8) 

The transformed equations (7) can, in principle, be solved for the 
phases pt(t) for t ^ 0. The desired stability information can be obtained 
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directly from these equations. We shall, however, obtain this informa- 
tion in a somewhat indirect but more profitable way by defining an 
associated linear time-invariant system with N inputs and N outputs. 

Consider the same control arrangement for the TV interconnected 
oscillators. Instead of the former starting conditions, suppose the con- 
trol paths have always been connected and that each oscillator can be 
activated by an arbitrary frequency "input" as shown in Fig. 2. Then 
the actual frequency of the ith oscillator at time / is the sum of the 
basic frequency input Vi(t) and the correction component r,(i) leaving 
the filter. The phases p,(t) are considered the "outputs" of the linear 
system. When y,-(0 — for each i, the system is at rest and all outputs 
Pi(t) are identically zero. 

The importance of the associated linear system is that any desired 
starting conditions in the physical model can be treated by an equiva- 
lent set of inputs to the linear system. To clarify this, note that the 
system of Fig. 2 is characterized by the equations 

pi(t) = ki(t) * E Oii\Pi(t ~ T tj ) - Pi(t)} + Vi (t), - co < t< co. (9) 



Formally taking the exponential (two-sided Laplace) transform of (9) 
we obtain 



sP< = Hi E da Pj - HtPi + Vi , 



(10) 



where Pi(s) and F,(s) are, respectively, the exponential transforms of 
Pi(t) and Vi(l). Equation (10) implicitly characterize the associated 
linear system whose inputs are v { (t) and outputs p,(i) as long as t>t(t) 
has an exponential transform. Comparing (7) and (10) we see that the 
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Fig. 2 — Model of station i of the associated linear system. 
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phase responses of Ihe physical model with the corresponding starting 
conditions will be the outputs of the associated linear system for t ^ 
if we select the inputs to be 



VM -i/< + p«(0) + QM. 

In the time domain, these inputs are 

Vi(t) = fm(t) + Pi(0)S(t) + qi (t), 



(11) 
(12) 



where qt(t) is the response of the filter Hi(s) to a time-limited input 
which begins at time t = —max r, 7 and ending at t = 0, 8(1) is the unit 
impulse function and u(t) is the unit step function. From (8) and the 
final value theorem it follows that Qi(t) —*■ as t -* °o . It is important 
to note that the phase responses to the inputs (12) will be the same 
as the phase responses of the physical model only for t ^ 0. For t < 
the responses of the associated linear system do not correspond to the 
physical model. 

Equation (10) may be expressed in the form 



Pi = &(«) E duPj + ( — i 



v if 



where 



his) = 



H f (s) 



S + Hi(8) ' 



(13) 



(14) 



The simplified model of the linear system, corresponding to (13), is 
shown in Fig. 3 where j8,-(s) is the transfer function of the feedback 
configuration as shown. Thus, the operation of the ith station is to 
average the incoming phases, apply the average to the filter /3»(s), 
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Fig. 3. — Simplified model of station i for the associated linear system. 
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and finally add a filtered input component to produce the phase response 
Pi(t). We shall see in Section IV that the condition for stability of the 
system is simply that the filter fii(s) have a gain less than unity for 
sinusoidal inputs. 

Equations (13) for i = 1,2, ■ ■ ■ , N can be formally solved for the 
phase responses with the help of matrix notation. Let B(s) be the 
.V X N matrix whose yth component is 

biiis) = fas) dij 

and let C(s) be the diagonal N X N matrix whose nth element is 

1 



cu(s) = 



s + H { (s) ■ 



Note that for s = we have c,,(0) = 1/X, and B(0) = A, where A, 
is the dc gain of //,-(s) and .4 is the averaging matrix, both defined 
earlier. Let P(s) and V(s) be the N component column matrices whose 
z'th elements are, respectively, P,(s) and Vj(s). Then (13) becomes 

[I - B]P = CV (15) 

or 

P = KV, (16) 

where 

K(s) = [I - B(8)}r l C(s) (17) 

is the matrix transfer function of the linear system. Thus, each element 
K,j(s) of K(s) is the scalar transfer function relating the output pi(t) 
to the input Uj(t) when all other inputs are zero. In Section IV we shall 
determine certain key properties of the singularities of K,j(s). In order 
to examine the behavior of / — B(s) in the neighborhood of s = 0, 
certain important properties of the averaging matrix A will be needed. 
In the next section these properties are derived. 

III. PROPERTIES OF THE AVERAGING MATRIX 

As a result of (4), the averaging matrix .4 has row sums equal to 
one. From the requirement that the network be connected, certain 
restrictions are placed on which combinations of elements of A may be 
zero. These two characteristics of A imply certain essential properties 
of the matrix / — .1 where / is the identity matrix. 
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Theorem 1: If A is the averaging matrix of a connected network of N sta- 
tions then the matrix I — A has rank N — 1 . 

Proof: Since the equation 

(/ - A)x = (18) 

is satisfied by any vector x with all components equal, the matrix I — A 
is singular and must, therefore, have rank less than N. Suppose that its 
rank is less than N — 1. Then (18) has at least two nontrivial solutions 
that are linearly independent. Therefore, there exists a nontrivial solu- 
tion, u, whose components are not all equal. Now let w be the solution 
vector with each component equal to the negative of the smallest 
component of u. Then y = u + w is a nontrivial solution of (18) with 
all components non-negative and at least one component equal to 
zero. Let G, = \i\ , it , • ■ • , n-\ be the set of indices for which y t = 
and 03 = [4+i , t*+i , • • • , is] the set of indices for which y t is positive. 
Since y satisfies (18) we have 

V 

Vi ~ YtOifilj = i = 1,2,---,N 

and so 

2 o-aVi = 0, for i e a. 

But this is only possible if 

a ,-j = 0, i e CL and j e (B, 

which implies that there is no transmission path from any station with 
identifying index in (B to any station with identifying index in Q. Con- 
sequently, the network of A r stations is not connected, which is a con- 
tradiction. Therefore, I — A must have rank N — 1. 

Theorem 2: If A is the averaging matrix of a connected network, the co- 
factors of all the elements in any given row of I — A are equal and posi- 
tive. Specifically, if M u is the cof actor of the ijth element of I — A, then 

Ma = M ik > 

fori,j,k= l,2,---,N. 

Proof: Since I — A has rank A r — 1, the solutions of 

(/ - A)y = 
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satisfy 



Vi = Ma 
y k M ik ' 



i = 1&>--,N. 



But the only solutions y are those with all components equal. There- 
fore, 

Ma = M ik all i,j,k. (19) 

Let R(c) = I — i'A and let M,j(p) be the cofactor of the ijth element 
of R(p). For S p < 1, each principal minor of R(c) is the determinant 
of a diagonally dominated matrix (see Appendix), so that 

Mu(£) ^ 0|c<l. 

Since il/,-,(0) = 1, it follows by continuity that M tl (l) ^ 0. Hence, 
from (19) 

Mij = M ik ^ all i,j,k. (20) 

Now (/ — .4)', where the prime denotes the transpose, must also 
have rank AT — 1 . Thus, solutions of 

(I - A)' z = (21) 

satisfy 

-'-£* i=l,2,---,N. (22) 

z k Mki 

Equal ions (20) and (22) imply that the nonzero components of z must 
have the same sign. Suppose a solution z of (21) has at least one com- 
ponent zero and nonzero components positive. Then the same argument 
used in Theorem 1 leads to the conclusion that the network is discon- 
nected, which is a contradiction. Therefore, there is a solution z with all 
components positive and consequently (22) implies that all cofactors 
M/j are positive, which completes the proof. 

IV. ANALYSIS 

With the help of the preceding results, we are now in a position to 
prove stability and determine the expression for final frequency. These 
results will be obtained under the assumption that /3,(s), for each i, 
satisfies the condition 

|/3.-0'«)l < 1, « * 0. (23) 
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In Appendix B we show, with the help of the Nyquist criterion, that 
condition (23) implies the stronger condition 

| j8,-(s) | < 1 for s in <R, (24) 

where (R is the right half and imaginary axis of the s plane excluding the 
point s = 0. 

In Section II we saw that the associated linear system is characterized 
by the matrix transfer function K(s) given by 

K(a) = [I - B(s)]-* C(s). (25) 

Now, since /3, = Hi/ {a + H t ), it follows from (24) that 

Cii (s) = 1/(8 + Hi) 

has no singularities in (R. Furthermore, under condition (24) the matrix 
/ — B(s) is diagonally dominated (see Appendix A) for all s in (R. 
Thus, the determinant | / — B(s) | is nonzero for all s in (R, and so we 
conclude that each component transfer function Kij(s) is analytic in 
the region (R. 

At s = 0, the matrix / — B(s) reduces to I — A which is singular 
according to Theorem 1. Thus, the determinant \I — B(s)\ has a 
zero at s = 0. To show that it is only a simple zero we find an asymptotic 
expression* for the determinant in the neighborhood of s = 0. In the 
matrix / — B(s), we replace the elements b</(«) by their asymptotic 
expressions 



s) ~ Oi 



1 " S T<j + t) S J ' 







where we have used the relations exp (— ar) **** 1 — st and H,/(s 4- Hi) 
r+j 1 — s/Xi . Without changing the value of the determinant, we may 
replace the first column by the sum of all the columns. The first column 
then becomes 



(-0 



where 

1-1 

is an average of the transmission delays of links arriving at the »th 



The technique for finding the asymptotic expression is due to A. J. Goldstein. 
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station. Now we expand the determinant about the first column and 
obtain 



I - B(s) |~sE (u + ~) M n , s -> 0, 



(27) 



where My is the cof actor of the ijth element of / — B(0), as defined 
in Section III. 

Since Ma is positive (from Theorem 2), it follows from (27) that 
| I — B(s)\ has only a simple zero at s = 0. But c,,(0) = l/\- is finite, 
so that from (25) we conclude that each Kij(s) has a simple pole at 
s = 0. Using (25) and the asymptotic expression (27) it follows that 

KM ~ yj/s, s -> 0, (28) 

where 

Mji/Xj 

7i = 



~Z(n + QM n 



Note that y, is positive and independent of i because My,- = M,i > 0, 
according to Theorem 2. Thus, letting 

we have 

1 * = 1 + ZndA, (30) 

with < dj < 1 and X) rf y = 1- 

We have, therefore, shown that each transfer function Ka(s) is 
analytic in the right half and on the imaginary axis of the s plane ex- 
cept at s = where it has a simple pole with positive residue inde- 
pendent of i. The impulse response kij(t), associated with K i} -(s), will, 
therefore, consist of exponentially decaying sinusoids and a step function 
of height 7/ . 

To determine the stability of the original model under the particular 
starting conditions, we examine the asymptotic behavior of the phase 
responses of the associated linear system when subjected to the inputs 
given by (11). From (16), we have 



PM = E km 



£/* + Pi(0) +Qi(a)l. (31) 



1700 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1966 

From the known properties of Kij(s), it follows that the phase response 
Pi(t) for t > will be the sum of terms decaying exponentially to zero 
plus a term of the form ft + ??,• where/ and jj, are obtained by the residue 
theorem according to 



and 



/ = lim f?Pi(s) 
*-»o 



* = lim£[« f /\00]. 
j-o as 



The final frequency / of the *th oscillator is, therefore, given by 

/ = £ Ti/i (32) 

2=1 

which is independent of i. Thus, we have proved stability of the system 
since the frequency of each oscillator has been shown to asymptotically 
approach the common frequency / and the phase differences clearly 
approach finite values. From (30) the expression for the final frequency 
can be written as 

£ djfj 

f = rrr ii (33) 

i 

which, with the help of (29) , shows the dependence of / on the delays 
nj and the dc gains A,- . 

V. REMARKS ON THE FINAL FREQUENCY 

From the results of the preceding section it is clear that the final 
frequency can be below even the lowest oscillator free-running fre- 
quency. In fact, it is evident from (33) that the final frequency is a 
monotonically decreasing frunction of the system gain-delay products. 
Thus, the controls may bring the system to a frequency outside its 
practical operating range. 

The final frequency reduction is a consequence of the fact that the 
frequency control of each station varies directly with the differences of 
total phase. The interstation delays introduce phase lags which drive 
down the frequency of each station. This point is made somewhat 
clearer by considering a system in which all the oscillators have the 
same frequency / and the same initial phase. When the controls are 
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applied at t = an average phase "error" (— fn) is applied to the con- 
trol path of each oscillator i. This "error" causes a simultaneous re- 
duction in the oscillator frequencies from which the system never 
completely recovers. 

As a conceptual solution to this difficulty, suppose the system of Fig. 
1 is modified so that the local phase at the output of the integrator is 
passed through a delay line before being compared with the incoming 
phases. This local delay at station i is chosen equal to r t - , the average 
delay of links terminating at station i, as defined in (26). 

In the previous model, the error signal was determined by a com- 
parison of the local phase at the present time with the remote phases 
of earlier times. In this modified system, however, the comparison is 
made between phases which on the average occur at the same time. 
Thus, the undesired component of the error signal due to interstation 
delays is eliminated. 

Using an argument which parallels the development of Sections II, 
III, and IV, it may readily be shown that the final frequency for the 
modified system is given by 

f=T,djfj. (34) 

i 

In contrast with the original system, it is evident that the final fre- 
quency of the modified system is always an average of the free-running 
frequencies. 

The BeneS formulation (Fig. 1) may be viewed as a simplified ab- 
straction of the more complex practical systems that have been pro- 
posed. 2 - 5 Both the Benes formulation and the modified system contain 
a total phase comparator which is an impractical element. Karnaugh 8 
has shown that an important linear subclass of the more realistic class of 
systems he has proposed obey equations of the same form (6) as in the 
Benes model. This more realistic formulation also fits the linear system 
model with modified frequency "inputs" that depend in a different 
manner on the initial conditions. It is, therefore, subject to the stability 
condition (23). 

Moreover, it has a different final frequency which approaches an 
average of the free-running frequencies as the interstation delays become 
large. 

In short, the Benes formulation was sufficient to provide the im- 
portant stability criterion, but neglected factors affecting the final 
frequency. The linear system model developed here is general enough 
to be applicable to both the Benes formulation and a linear subclass of 
the more realistic Karnaugh formulation. 
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"VT. SENSITIVITY AND RELATED QUESTIONS 

Suppose that the system has been operating in a synchronized steady- 
state condition for a long time, and at some instant, say t = t o , a sudden 
change is made in one or more parameters of the system. The subse- 
quent phase responses are determined by considering the new associated 
linear system subjected to suitable inputs equivalent to the pertinent 
starting conditions. These inputs will have exactly the same form 
as (11) but now the term Q,(s) will be evaluated using the past history 
of the phases given by 

Pi(t) = f s t + m t<U, 

where /, is the synchronous frequency prior to the parameter change. 
If the stability condition (23) is satisfied for the new system and if the 
parameter change does not reduce any X,- or a i} - to zero, the new system 
will also be stable. Consequently, after t = to the frequencies of the 
oscillators will asymptotically resynchronize to the new final frequency 
determined by (33) or (34) using the changed parameter values. From 
these arguments we can also deduce that the effect of a slowly time- 
varying parameter on the system operation is to cause a corresponding 
slowly varying synchronous frequency. By "slow" time variations we 
mean that the time for a noticeable change in a parameter value to 
occur is much longer than the time constants associated with the tran- 
sient response of each i£#(s). 

By similar arguments, it is easily seen that failure of a transmission 
link will lead to resynchronization if the remaining network is still 
connected. Also, in the case of oscillator failure, the remaining N — 1 
oscillators will resynchronize to a new frequency if the resulting network 
of N — 1 stations is still connected after removal of all transmission 
links entering or leaving the inoperative station. In each case, the final 
frequency can be computed from (33) or (34) using the appropriate 
parameter values. To prove these results, the nonzero averaging weights 
?an be rescaled so that A has row sums unity; the filter gains X,- are 
assumed to be correspondingly rescaled. The characterizing equations 
for the new system then has the required form and so resynchronization 
vill occur. 

The effect of independent jitter noise on the frequency of each oscil- 
lator may be considered by including a noise term n,-(0 in each "input" 
Vi(t). By superposition, the effect of noise can be considered separately. 
Thus, each phase response pi(t) will consist of the response in the 
absence of noise plus a noise component whose power density spectrum 
is 
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y=l 

where &(«) is the power density spectrum of n,(t). Consequently, if 
the input noise jitter has zero mean and finite variance the output 
noise components will also have zero mean and finite variance. We 
conclude, therefore, that in the presence of noise jitter each oscillator 
will asymptotically have a common frequency with a random perturba- 
tion. The perturbations will be correlated but, in general, will not be 
identical. Furthermore, small jitter noise implies proportionately small 
perturbations. 

VII. CONCLUDING REMARKS 

We have seen that the transfer function approach has permitted a 
simple treatment of a rather complicated control system. Further 
studies regarding transient response or bounds on the size of perturba- 
tions due to jitter noise can be made for particular topological configura- 
tions by determining more information about the transfer functions 
Kij(s) with the help of (25). The linear system approach together with 
the added generality of having different filters at each station has 
made it possible to consider the effect of parameter changes or oscillator 
failure on the behavior of the system. 
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APPENDIX A 

A square matrix A is said to be diagonally dominated when for each 
row the sum of the magnitudes of the off -diagonal elements is less than 
the magnitude of the diagonal element, i.e., 

I an I > £ I By I » eacn i- 

Theorem: If A is diagonally dominated it is nonsingular. 

Proof: Suppose the contrary. Then there exists a nontrivial solution 
[Xi] satisfying 

S dijXj = each i. 



1704 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1960 

Let r be one of the indices for which | .r,- | is a maximum. Then 

a rT X T = — 2 a rk%k 

k*T 

so that 

I «rr | | X r | ^ 2 | O-rk \ \ X k | ^ £ | d rk \ \ X T \ 
k^r k^r 

which is a contradiction. Hence, the theorem is proved. 

APPENDIX B 

On the Boundedness of /3,(s) 

Theorem: If @i(s) is bounded by unity on the joi axis then it is bounded 
by unity in the entire right-half plane. 



Proof: Since 



PM ~ l+Hds)/*' 



the condition | /3,(jo>) | < 1 is equivalent to 

\A\ < | 1 + A exp (i<p)\, (35) 

where 

A exp (up) = Hi(ju)/ju. 

But (35) is equivalent to 

A cos ip > — \ 

so that the locus of H(j<a)/j<a, as o> increases from — oc to », cannot 
encircle the point — 1. Hence, by the Nyquist stability criterion, /3,(s) 
is analytic in the right -half plane. Furthermore, /3,(=o) = 0. Thus, 
it follows that |j?*(s)| < 1 in the right -half plane according to the 
rmiximum modulus theorem. 
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